Ordering in spatially anisotropic triangular antiferromagnets 
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We investigate the phase diagram of the anisotropic spin-1/2 triangular lattice antiferromagnet, 
with interchain diagonal exchange J' much weaker than the intrachain exchange J. We find that 
fluctuations lead to a competition between (commensurate) collinear antiferromagnetic and (zig- 
zag) dimer orders. Both states differ in symmetry from the spiral order known to occur for larger 
J', and are therefore separated by quantum phase transitions from it. The zero-field collinear 
antiferromagnet is succeeded in a magnetic field by magnetically-ordered spin-density-wave and 
cone phases, before reaching the fully polarized state. Implications for the anisotropic triangular 
magnet CS2CUCI4 are discussed. 
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The nearest- neighbor spin-1/2 Heisenberg antiferro- 
magnet on the two dimensional triangular (hexagonal) 
lattice has long been a prototypical theoretical model 
for frustrated quantum magnetism |IJ. A simple and 
relevant generalization is the nearest-neighbor spatially 
anisotropic triangular antiferromagnet, defined by 



H 



(1) 



Here i,j are sites on the triangular lattice, Jij — J, J' 
along horizontal and diagonal links, respectively, and 
zero otherwise (see Fig. [TJ). The quasi-one-dimensional 
(ID) inorganic salt, Cs 2 CuCl 4 [2], with J'/J = 0.34, pro- 
vides an experimental application. Considerable theoret- 
ical interest in CS2CUCI4 was spurred by the observation 
of a strong inelastic continuum in neutron scattering^, 
prompting many interesting two-dimensional spin-liquid- 
based interpretations 0, [1, E|- There are also very 
detailed studies of the magnetization process and as- 
sociated temperature (T) vs magnetic field (h) phase 
diagram^, 0]. Existing quantitative comparisons with 
experiments are based on semiclassical spin-wave calcu- 
lations 0,0, This well-studied scheme cannot, how- 
ever, be applied to the very interesting (and, we argue be- 
low, experimentally relevant) Id limit of J' <C J, where 
the system is best described as a collection of weakly 
coupled quantum-critical spin-1/2 chains. 

In this Letter, we report a systematic analysis of this 
important limit, as a function of magnetic field, includ- 
ing the interesting case of zero field. This is accom- 
plished using a Renormalization Group (RG) analysis 
of the perturbations, represented by the interchain ex- 
change J', to the exactly solved problem of decoupled 
Heisenberg chains. As the RG progresses, new inter- 
chain couplings, consistent with symmetries of the lattice 
model ([TJ, are generated. According to standard RG ar- 
guments, the low energy physics is controlled by the cou- 
plings which renormalize to dimensionless values of order 
one first. We find that in zero field all dominant inter- 




FIG. 1: (color online) Left: The triangular lattice, showing 
(J, J', J2) exchange interactions, and integer coordinates x,y. 
Black (red or gray) arrows indicate spin order in the CAF 
phase on odd (even) chains (J2 = section). Solid ovals 
show the strong bonds in the ordered zig-zag dimer state 
(J2 = 0.241 section). Right: suggested zero-field phase di- 
agram (solid points indicate position of phase boundaries on 
the axes). 



chain couplings are generated by fluctuations, and missed 
in a naive analysis. The fluctuation-generated relevant 
interactions foster two competing orders: collinear an- 
tiferromagnetism (CAF) and spontaneous dimerization. 
For the simple Heisenberg model, this competition is re- 
solved by a marginal intra-chain interaction, in favor of 
the CAF state. A more general picture, illustrated in 
Fig. [TJ emerges in the presence of an additional second- 
neighbor exchange J2 along the chain axis, exposing the 
competing (four-fold degenerate) spontaneously dimcr- 
ized zig-zag state. In a magnetic field, we recover the 
main low energy features of the CS2CUCI4 experiments. 

We begin with the zero field case. The low energy prop- 
erties of each decoupled Heisenberg chain are described 
by a critical Wess-Zumino-Novikov-Witten (WZNW) 
SU(2)i theory, with central charge c = 1. The low energy 
theory is connected to the lattice model by the continuum 
limits of the spin and exchange operators: 



My{ X ) + (-l) X N y ( X ), (2) 

(-lfe y {x). (3) 
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The dominant fluctuations of a single chain are evident 
above in the operators Ny(x) and e y (x), which both have 
scaling dimension 1/2 and stronger correlations than any 
other operators. We have also introduced the the con- 
served uniform (spin) magnetization density M y (scaling 
dimension 1). In terms of the usual right (R) and left (L) 
moving chiral spin currents, M y = Jr >v + Jh, y - Here we 
use the integer coordinates x, y indicated in Fig[TJ 

Let us now consider the operators which perturb the 
WZNW theory. This includes a marginal "backscatter- 
ing" term already for a single Hciscnbcrg chain, and, in 
addition, terms arising from inter-chain coupling. In un- 
frustrated systems of coupled chains, e.g. a simple rect- 
angular lattice, a naive treatment of these latter interac- 
tions, that consists in replacing the spin operators in the 
inter-chain exchange J' using Eq. Q , is sufficient ll| . In 
our case this produces the following perturbation: 



H' = J dx {jb S Jn,y ■ Jl, 



(4) 



+-f M M v ■ My+1 + ItwNy ■ d x N y+1 }, 



with 7bs ~ O(J), 7 M = 2J', 7t w = J'. The effect of 
frustration is evident in Eq.(j4]) by the absence of direct 
N y ■ N y+ i coupling between chains. Instead, one obtains 
only the small residual "twist" term 7 tw , whose scaling di- 
mension is increased by 1 by the spatial derivative. This 
marginally relevant term leads on its own (for instance in 
a two-chain zig-zag ladder) to a very weak incommensu- 
rate spiral spin order with an exponentially small gap.[l2j| 
However, the naive interactions in Eq.([4|) are incom- 
plete, due to the dynamical generation of strongly rele- 
vant interchain interactions. From simple power count- 
ing, the most relevant interactions possible should involve 
products of the N y or e y operators on different chains. 
Between nearest-neighbor chains, however, such interac- 
tions are forbidden by symmetry under reflection in a 
plane perpendicular to the chains. However, there is no 
symmetry prohibiting relevant couplings between second 
neighbor chains: 



H" = ^2 fdx{g N Ny-Ny + 2 + g e e y £y + 2}. 



(5) 



The terms g N , g £ have scaling dimension 1, and if present, 
are strongly relevant. Being symmetry allowed, one ex- 
pects them to be generated by fluctuations (in the RG 



sense) on short scales, similarly to Refs. 

A convenient method to derive (O from (j4]) is to for- 
mally perturbatively integrate out, say, odd chains from 
the action corresponding to (|4]). While the resulting 
even-chain action is naturally non-local, it is conveniently 
treatable by the RG. At first order in this action, the 
RG generates a local derivative term oc d x N y ■ d x N y+ 2- 
At second order, the RG induces the terms in (O 
by more standard calculations very similar to those in 



Refs. 13j, This happens already during the initial 
stage of the RG (0 < I < £ ~ 0(1)). Because these 
terms then grow rapidly with increasing i, further gen- 
eration of these terms can be neglected. Therefore the 
effect of fluctuations is to determine "initial conditions" : 

id \ ^ , . 3^4g7jv/7tw 3^4o(J') 4 (ei \ 

= — 2 g»{io) - -^r- = (6) 

using the values of 7 tw , j M and velocity v = nJ/2. Here 
Aq is a normalization factor for the N field [l5| . Despite 
the smallness of g e {to),g N {lo), it is clear that they will 
dominate the low energy physics, since the scale required 
for the remaining marginal interactions to grow large, 
L 



marg 



exp[c( J/ J')] is exponentially long. 



Naively, the competition between g e , g N is keen, since 
the generated values in Eq.© are of the same order. 
However, the outcome is definitely determined by the 
marginal intra-chain backscattering interaction 7b s . This 
can be seen by considering the coupled RG flow of these 
three interactions (we neglect terms of 0(J' 6 ) or smaller): 



dtjhs = 7b s 



1 3 

In = g N - ^7bs<?iv, dig e =gs + ^7bs5e,(7) 



where gx = gx/{2nv) (X = bs,N,e). Here % s (£ = 
0) = -r w -0.23 < 0,[lf| so the backscattering is 
marginally irrelevant and flows logarithmically to zero: 
7bs(^) = — r/(l + T£). For an isolated Heisenberg chain, 
7bs only weakly modifies the low energy behavior by loga- 
rithmic corrections [171 | . Because 7b s enters the RG equa- 
tions for g e multiplicatively, however, it becomes cru- 
cial in the coupled-chains problem. Indeed, from 1(7} the 
relevant couplings grow according to 



9n{£) = „l/4 i- 



(8) 



where E e = {1+T£)/(1+T£ ). With the small initial val- 
ues in Eq. (J6]) , the relevant couplings become of order one 
at the scale I* - 41n(J/J') > 1. Then since Yt > 1 we 
have |"^f, ^ — 3> 1. Thus backscattering parametri- 
cally enhances the antiferromagnetic interaction between 
second neighbor chains over the competing dimerization 
one. As a result, subsystems of even and odd chains or- 
der into Neel patterns independently of each other on the 
scale £*, when g N (£*) ~ 1. The remaining coupling j M 
between these two rectangular Neel sublattices remains 
small and can still be treated perturbatively. By a stan- 
dard quasi-classical spin wave calculation, appropriate 
for the magnetically ordered Neel sublattices, one ob- 
serves the familiar phenomenon of order-from-disorder, 
which locks the two sublattices into the collinear CAF 
phase shown in Fig. [TJ For values of J' w J, it is 
known that the ground state is a rather different spiral, 
smoothly deformed from the famous 120° order obtain- 
ing for J' = J [3]. These two phases differ in symmetry, 
and must be separated by at least one quantum phase 
transition (which can be 1 st order) when J' /J ~ O(l). 
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Now consider the effect of an additional second- 
neighbor antiferromagnetic exchange along the chains, 
J 2 . This is well-known to decrease the initial value of 
the backscattering, T. In particular, T — * at the criti- 
cal value J" 1 * = 0.241 [16j, which separates the critical 
(Luttinger) phase of the spin chain from a spontaneously 
dimerized one. With non-zero J', sufficiently close to this 
critical point one finds T£* ~ 1, and the backscattering- 
induced enhancement of g N is eliminated. From Eq.([6]), 
the larger initial value of \g e \ induces instead a two- 
dimensional dimerization instability! Moreover, compar- 
ing the RG scales of the relevant, £*, and backscattering, 
&bs — 1/r, terms in the vicinity of the J£ rlt critical point, 
we deduce the phase boundary between the collinear and 
dimerphases: J'/ J| C -d ~ exp[-0.11J/( J 2 crit - J 2 )]. Sym- 
metry again requires at least one additional phase bound- 
ary (dashed in Fig. [T]) between the spiral and dimerized 
phases. The nature of the transitions in Figfl]is beyond 
the scope of this paper. 

Let us turn now to the situation in non-zero field. 
It is convenient to work at fixed magnetization, M — 
~ ~k^2y Jd x My(x), rather than fixed field. 
The Id Heisenberg model retains c = 1 free-boson char- 
acter for any < M < 1/2, and can be viewed as an 
"easy-plane" deformation of the WZNW model. All scal- 
ing dimensions can be expressed in terms of the "boson 
radius" R, which is a known function of M [19]. For 
< M < 1/2, the transverse XY components of the 
staggered magnetization field = N x ±iN y strengthen 
their correlations, and have scaling dimension A xy = ttR 2 
decreasing from 1/2 at M = toward 1/4 at M = 1/2 - . 
This is in accordance with the semiclassical canted XY 
antiferromagnetic order in a field. Conversely, the stag- 
gered magnetization along the field direction shifts to the 
incommensurate wavevectors tt ± 25, 

N*(x) S z y + e a5x + S*-e~ mx , (9) 

with 5 = ttM. The correlations at this wavevector pro- 
gressively weaken as M is increased, such that S z± has 
scaling dimension A z = (AttR 2 )^ 1 increasing from 1/2 
toward 1 as M — > 1/2 _ . The dimerization operator e y is 
expressible as a different linear combination of S z± and 
is not an independent degree of freedom for M ^ 0. Fi- 
nally, the zero field magnetization operator develops an 
expectation value: My(x) — > M + 5My(x), and the fluc- 
tuation piece 8My(x) has scaling dimension 1 as before. 

These changes in the correlations have several effects. 
The shift of the longitudinal spin correlations to the in- 
commensurate wavevector 7r ± 25 removes the frustration 
of the longitudinal piece of the inter-chain J' spin interac- 
tion. The Z-components of neighboring spins on a chain 
are no longer antiparallel, and the effective exchange field 
exerted by them upon the shared spin on the adjacent 



chain does not cancel. This leads to a relevant term 

#sdw = E [ dx Tsdw^+S^x + h.c, (10) 
y J 

with 7 s dw = (1 — e 2lS )J' non-zero for M > 0. For small 
magnetization, this term is almost as relevant as the 
fluctuation-generated terms in Eq.([5]), but its magnitude 
is O(J'S) and so is significantly larger for 5 > (J'/ J) 3 . 
For field above h% « 3 J 4g/(87r)(j v / o f) 3 ' l sat ^ becomes the 
dominant instability channel, and stabilizes a longitudi- 
nal "spin density wave" (SDW) with the incommensu- 
rate wavevector k x — 7r ± 28 and all spin expectation 
values aligned along the field. However, as M is further 
increased, its scaling dimension A SDW = (2TTR 2 )- 1 be- 
comes less relevant, approaching 2 at the saturation. 
Simultaneously, the XY piece of the twist term, 

HL = E ! dx ( N v d * N v+i + h - c ' ( U ) 
v J 

becomes increasingly more relevant with M (as com- 
pared to marginal in zero field). Its scaling dimension, 
A tw = 2ttR 2 + 1, decreases steadily from 2 toward 3/2 as 
M — > l/2~. The transition between two ordered states 
takes place when the two scaling dimensions become ap- 
proximately equal at irR 2 = (\/5 — l)/4 ~ 0.309, which 
corresponds to M « 0.3 (see Figures 8 and 9 in Ref . fioh . 
Thus for M > 0.3, the twist term dominates and in- 
duces a "cone" state, in which the XY components of 
the spins spiral at a wavevector close to but not equal 
to k x = 7T. One may wonder whether the fluctuation- 
generated terms in ([S]), particularly the transverse piece 
of g N , can intervene between the two described states. 
However, because it is only generated at fourth order in 
J', we find that it is always subdominant to either 7 s d w 
or 7t w , at all fields. The resulting phase diagram of J—J' 
model in magnetic field is sketched in Figure [2] 

We now apply this understanding to CS2CUCI4 . We 
first address an objection to the quasi-ld approach, 
raised in Ref. 0, that the magnetization curve M{h) of 
CS2CUCI4 differs significantly from that of the Id Heisen- 
berg chain. In fact, very good agreement can be obtained 
simply by assuming that spins are completely uncorre- 
cted between different chains. By calculating the expec- 
tation value in a direct product state on different chains, 
the ground state energy suffers a correction obtained by 
simply replacing J'Si-Sj — > J'M 2 for each J'-link. Using 
h = —dE/dM and the known thermodynamics of decou- 
pled Heisenberg chains, one obtains the curve in Fig. [2] 
Thus the observed M(h) curve in fact indicates the weak- 
ness of inter-chain spin correlations in CS2CUCI4. How- 
ever, by its definition, M{h) depends only upon (Si ■ Sj) 
for nearest-neighbor This explains how both spin- 
wave theory[3, IH and our quasi-ld approach, which dif- 
fer drastically in longer-range correlations, can reproduce 
the experimental M(h) curve. 
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FIG. 2: (color online) Left panel: Comparison of experimen- 
tal (solid black line - obtained by scanning Fig. 2a, B || a, 
of Ref.Q) and theoretical (dashed red line) magnetization 
curves for CS2CUCI4, J'/ J = 0.34. The experimental (and 
theoretical) M(h) is nearly direction-independent apart from 
(/-factor rescaling. Right panel: T — phase diagram of J — J' 
model (top line) and CS2CUCI4 (h || a and h _L a lines). h sa .t 
is the field at which M = 1/2 is reached. 



Long-range correlations are crucial for understand- 
ing the low-temperature phase diagram of CS2CUCI4. 
An important (and unfortunate!) complication of this 
material is a significant Dzyaloshinskii-Moriya (DM) 
coupling along the interchain links. Although nom- 
inally small, D w 0.05 J f» 0.16 J', it plays a key 
role because it is non-frustrated. Specifically, apply- 
ing Eq.([5]) gives the inter-chain DM coupling £?dm = 
J2 y J dx 2D {N^N^ +l - NyN* +1 ) , taking conventional 
crystalline axes a, b, c. In zero field, this term is as rel- 
evant as the fluctuation-induced interactions in ([6]), but 
is much larger since D/J ^> (J'/ J) 4 . Thus it dominates 
and drives spiral magnetic order, as seen experimentally. 
Magnetic field along the a axis (h || a in Fig. [2j) strength- 
ens {b, c} components of N, and, hence, the DM-cone 
phase, which extends all way to /i sa t- 

When the field lies in the b — c plane (h _L a in Fig. \2§ , 
the situation is more interesting. In this case, the DM 
coupling involves both XY and Z components of the TV 
field. From ^ we see that the correlations of these two 
fields become incommensurate. This effectively nullifies 
the DM term once the incommensurability 25 > AqD/J. 
Thus most of the behavior in this field orientation can be 
understood from the simple J — J' model in a field (top 
line in Fig. [21 right panel) discussed above. The same 
conclusion was reached in Ref.@. The experimental phase 
diagram for h _L a in Figj2] of Ref. 7} , indeed shows evi- 
dence of the expected DM, SDW, and cone phases. More- 
over, in neutron measurements [2j, the ordering wavevec- 
tor is consistent with the 5 = ttM expected for the SDW 
state for H = 1 — 2T, and evidences the cone state near 



At fields intermediate between these two limits, how- 
ever, experiment in addition observes one or two (depen- 
dent upon field orientation in the b — c plane) commen- 
surate ordered states, [7|, [20( denoted as SDW(CAF) in 



Fig. [2l An explanation within the J — J' model was pro- 
posed in Ref. |20|, based upon an extrapolation of the spin- 
flip expansion about the fully-polarized state. Some par- 
tial confirmation of this notion is found in our approach: 
the symmetry of the observed order is consistent with 
what would be induced by the fluctuation-induced g N < 
interaction in a field. [2l[ However our calculations above 
indicate that the magnitude of g N is much too low to 
explain the experimental ordering temperature T c (/i)[3| 
(note there are no logarithmic enhancements analogous 
to those in Eqs.JSJ) in large fields). We instead suggest 
that a relatively small direct antiferromagnetic exchange 
interaction J' 2 between spins on second-neighbor chains 
is a more probable explanation. We find only J 2 w 0.06 J 
is needed to explain the flat T c (h) curve 
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